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ASYMPTOTIC BEHAVIOURS OF SOLUTIONS FOR 
FINITE DIFFERENCE ANALOGUE OF THE 
CHIPOT-WEISSLER EQUATION 

by 

Houda Hani and Moez Khenissi 


Abstract. — This paper deals with nonlinear parabolic equation for which a local solution 
in time exists and then blows up in a finite time. We consider the Chipot-Weissler equation: 

2p 

Ut = Uxx + uP - , X G (-1,1); t > 0, p > 1 and 1 < g < ——. 

p + I 

We study the numerical approximation, we show that the numerical solution converges to 
the continuous one under some restriction on the initial data and the parameters p and q. 
Moreover, we study the numerical blow up sets and we show that although the convergence 
of the numerical solution is guaranteed, the numerical blow up sets are sometimes different 
from that of the PDE. 
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1. Introduction 

In this paper, we consider the nonlinear parabolic problem 

{ Ut = Uxx + UP - \Ux\'^ , X e (-1, 1), t > 0, 

u(±l,t) = 0, t>0, (1) 


ICey words and phrases. — Chipot-Weissler equation, blow up, finite difference scheme, numerical 
blow up set, asymptotic behaviours, numerical convergence. 
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Here p > 1, 1 < q < 


2p 

p + l 


and Mo is a positive function which is compatible with the 


boundary condition. It is well known that for some initial data, this problem blows up 
in a hnite time. Problem Q was studied for the hrst time by Chipot and Weissler in 
[ 1 ], since then, the phenomenon of blow up for different problems has been the issue of 
intensive study, see for example HI, 01,0,171,0 and the references therein. There exists 
many theoretical studies on the question of the occurence of blow up, but from a numerical 
point of view, many interesting numerical questions for problem ([^ are not treated. 


We define the blow-up set for problem ([^ as: 

B{u) = {x G [—1,1]; 3 {xn,tn) —t {x,T*) such that u{xn,tn) —t +C )0 as n —>■ -foo} . 

It is proved in [2] that the solution of ([^ blows up only at the central point, that is: 

3 T* < -foo such that lim u(t, 0) = -|-cx) but lim u(t, x) < oo when x ^ 0. 

t^T* t^T* 

In |5], we have conctructed a finite difference scheme whose solution satishes the same 
properties as the exact solution and moreover, we have proved that its solution blows up 
in a hnite time. In this paper and for the same scheme, we show the convergence of the 
numerical solution to the continuous one under some restrictions on p and g, and we study 
the asymptotic behaviour of the solution near its singularity. We prove that the numerical 
solution can blow up at more than one point, while a one point blow up is known to occur 
in the continuous problem. More precisely, we show that even if a difference solution blows 
up, its values remain bounded up to the moment of blow up except at the maximum point 
and its adjacent points, moreover, the number of blow up points depends, in a way, on the 
value of the parameter q. 


We recall the scheme studied in [Fj, for j = 1,..., and n > 0 we have 




= Uo{Xj), 
Un = m” 


U 


i+1 


- 


"-n 


+ - 


{2hr, 


yi+i - 


^ _ 7/^+1 I 

ni+1 S-ih 


-/Vn + 1 


= 0 . 


We denote by f/” := (mq, the numerical solution of (|^, and 

\\U^' 


( 2 ) 


= max \u. 


the L°° norm of U^. 

Here the notation w" is employed to denote the approximation of u{xj,t"‘) for Xj G [—1,1] 
and > 0. Also, we fix other notations as follow: 


1. r : size parameter for the variable time mesh r^. 
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2. h : size parameter for the variable space mesh /i„. 

3. t": n-th time step on t > 0 determined as: 

r to = 0 

J n—1 

I tn = + Tn-l = Y. '<'k] n>l. 

\ k=l 

4. Xj\ j-th net point on [—1,1] determined as: 

{ Xo = -1, 

Xj = Xj_i + hn, j > 1 and n > 0, 

XNr^ + l = 1 . 

We suppose that a spatial net point Xm coincides with the middle point x = 0. 

5. Tn '■ discrete time increment of n—th step determined by 

T„ = r min (l,||f/"||J+‘). 


Q. hn'- discrete space increment of n—th step determined by 




7. iV„ 

8. m - 


- -1 the number of subdivisions of the interval [—1,1]. 

Nn + 1 


As in [5], we suppose that the initial data uo satisfies the following conditions: 
(Al) uq is continuous, nonconstant and nonnegative in [—1,1], 

(A2) Uq is spatially symmetric about x = 0. 

(A3) Uq is strictly monotone increasing in [—1,0]. 

(A4) «o(-l) = «o(l) = 0. 

(A5) Uq is large in the sense that UmoIIoo 


This paper is organized as follows: In section 2, we state and prove the main results, 
that is, if p = 2 and q = 1 then the solution blows up at the maximum point and the points 

2(p — 1) 

around it, but remains bounded at all of the rest points, while if p > 2 and q < -, 

p 

then there is only a single point for the solution to blow up. In section 3, we prove 
the convergence of the numerical solution to the exact one. In section 4, we give an 
approximation of the blowing-up time. Finally, in section 5, we present some numerical 
simulations. 
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2. Main theorems 


In this section, we study the asymptotic behaviour of the difference solution near the 
maximal point Xm- 

Theorem 2.1. — Let be a solution of ([^, we suppose that h < —. For p = 2 and 

q = 1, we have 


lim _i = lim , i = +oo. 

n^+oo n^+cxD 

Proof — For j = m — 1, the equation of (|^ can be rewritten as 


(1 + 2An)'Um-\ 


(3) 




T-n. 


{2h, 


Pm - «m-2 




- u. 


n+1 I 
m—2 


Using positivity and monotony we get 


(1 + 2a„)<a > A„<+' + 


n+1 I n 


Tr, 




( 2 + 


We use that 

we obtain 


< - <-2 < 2< and <+■ - < 2<+'. 


, n+1 ^ ■''n „ n+1 ' 

Wm-i-i^2A„™ 1 + 2A„“™-^ 


A. 


Tn 


-u^ 


Furthermore from (|^ for j = m, we have 


n+1 _ ^-^n ^ ,n+l 


7 / — — 

™ 1+ 2A 


«+-i + 




which implies that 


1 + 2A, 




hn(l + 2 An) 

(1 + r+Oa-'), 


y;" 

^m) 


n+1 ^ ^rn 

^ - 1 + 2A, 


(4) 

(5) 

( 6 ) 

(7) 


Using (|^, ^ and 0 we get for p = 2 and q = 1 




A. 


-u. 


n 

(1 + 2A„)^ 1 + 2A„ ^ ^ /i„(l + 2A„) _1 + 2A 


2A 


^ „,n+l 


U. 


Then, 


1 + 


2'Tn\r 


ul+K > 


An n 1 7 

+ 3- 


«r-\ + Y^(l + ^n«") 

TnK,{l + TnU'ff) 


hn{l + 2\ny) ^-^-{l + 2KY 1 + 2A„ hPl + 2A,,)2 ’ 

which implies that 

n+1 \ ^nhn'U^ + hn(l + 2An)Mm-l ~ '^n'Wm(7 + 


'^m-1 — 


+ (1 + 2 An)^ + 2TnXr, 


( 8 ) 
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Since the solution blows up, then we have > 1, moreover 


Tn = — and hn = min(A/ 2 , h) = h. 


Then 


Hence, ([^ implies 




r 




T T 

hn Ar) — - and u, At, = 


hul 


h‘^{u 


n W ' 
'mJ 


U 


.n+l 


'm-1 — 


> 


r / 2r 
-—h h I 1 + —— 

h V h^u'l 


u. 


m—1 


■ r(l + r) 


h I 1 + 


2r 


+ 


2r' 




As we have 


then we get 


lim vh = +CX), 


n^+cxo 


r 


- + h liin + 

lim <t\ > ^ -. 


n^+oo 

If we assume that lim uh, ^ +cxd, let I = lim then we have 

n-)-+oo n-)-+oo ™ ^ 

, T , r(l + r) r, 1 


1 + r 


which is a contradiction because h < 

Therefore, we have 

lim m” _i = + 00 , 

n^+oo 

and using symmetry we get the result of Theorem 2.1. 

The next important result for this paper is mentioned in the next theorem: 

Theorem 2.2. — Let U"' he the solution of (1^, we suppose that p > 2 
2{p-l) 


p 

(a) If p = 2 and q = 1 then 

lim _2 < Too. 

n^+oo 

2(n — It 

(b) If p > 2 and q < - - then 

p 

lim < + 00 . 

n^+oo 


□ 


and 1 < q < 
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Proof. — Let prove (a): In (|^, if we take p = 2, q = 1 and j = m — 2, we get 


u, 


n+l „,n 


m-2 ^m-2 _ ^m-1 ‘^m-3 , (^,n \2 A,n+l i,n+l 

^ “ fpi \^m-2) 2h 


n+l _ 27,’^+1 -L 7,«+l „ 

“m—1 2 “m—3 , A.n 

72 ^ \^m-2} ) 


< 


"'n 


bnt < 0) then 




m—2 “m—2 ^ “m—1 ^m—2 


< 


ul+\ - 


Tr,, 


"■n 


+ (<-2)'. 


which implies that 

(1 + -^n)Mmt.^2 ^ + (^ + '^nKn- 2 ) Kn-2- 

In the other hand, in (|^ if we take j = m — 1, we get 

7,”+l _ 7,^1 7,"+l _ 2?/"+^ -7 ?/”+^ 

“m-l “m-l ^ “m-2 ^“m-1 “r 


< 


Tn 


/l2 

“n 


+ «-i) . 


bnt vfff }2 — u^\ < 0, then 


7,"’+l _ 7,n —7,’^+! -U 7,”+l o 

1 ^m—1 ^ ^m—1 ^ , /^^,n 


Tr,, 


h2 

“n 


which implies that 


and then 


“m-l — 

Next, if we recall ([^ for p = 2 we get 

n+l ^ 

^ 1 + 2A, 


(1 + < A„m"+^ + (1 + T-n^m-l) “m-D 


„+l , + (1 + T-nW” _l) <_l 


” 7,’"+! -U 

“m-1 


1 + A„ 

1 + TnU. 


1 + 2A 


n“m n 


Pntting 01 in ( [l0| we get 


7 ,n+l <- 

“m-1 — 

which implies that 


A. 


1 + A„ 


^■^n n+l I 1 + TnUj^ 

I + 2 A 7 , 1 + 2A„ 


{l+TnU^_^) 

+ ^^- -U 

1 + An 


m—1 ’ 


2A? 


^n+i / A,,(l + r„0 ^ l + r„<_i ^ „ 

_1 ^ -x \ /-. •« \ ^TY} I - •« 


(1 + A„)(l + 2A„)/ (1 + 2A„)(1 + A„) 1 + A 


m—li 


and then 


7 , n+l <- 
“m-1 — 


A,,(l + TnUlfjUl^ + (1 + 2A,7)(1 + i 


( 9 ) 


( 10 ) 


( 11 ) 


1 + 3A,i 


(12) 


























ASYMPTOTIC BEHAVIOURS 


7 


Now, putting (12) in (|^, we get 


7/”+^ < 

— 


1 + 


A^(l + TnU^)u^ + (1 + 2A^)(1 + TnUl^_{)U^_^ 
1 + 3A„ 


+ 


(1 + r„<_2) <_2 

1 + An 


(1 + TnUl^_^) ^ A^(l +rnO < + An(l + 2An)(l + Tn^.i)?/” _i 


1 + Ar 


(1 + An)(l + 3An) 


Then 

here we have put 

and 


u'^-2 < AnU^_2 + Bn, 


An — 


(l + TnU^2) 


1 + A^ 


Bn = 


An(l + ^nKn)u^ + An(l + 2An)(l + rn<_iX_i 
(1 + An)(l + 3An) 


Then the inequality (13) implies by iterations that 
Kn_2 < Aln-in"-2 + ^n-1 


^ An-lAn-2U^\ + An-lBn-2 + -Bn-1 


(13) 


n—1 


n—2 


n—1 


^n-1 


^n-1 


^ “m-2 ^ yBk Aij + Bn 

k =0 k =0 \ i = k-\-l J 

n n—2 n—1 

— '^m-2 Bk Ai + B„ 

k =0 k =0 i =0 

n n—2 n 

< u^_2 Ak + Bk Ak + Bn-i j A), 

k =0 k =0 k =0 

n n—1 n 

A;=0 k =0 k =0 

n n n 

< '“m-2 Y\ Ak + y^^BkY\ Ak 


\ k =0 


k =0 


k =0 k =0 


< <-2+e^'= n^‘ 


k =0 / k =0 

To ensure boundedness of u '^_2 we shall prove that 


Bn < +CXD and An < +cxd. 


n>0 


n>0 


To do this, we need the next lemma: 
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Lemma 2.3. 


We define the sequence an = 


u. 


m-l 


U. 


n 


1. For p = 2 and q = 1, we assume that supM!i_i > 717(1 + r), then (a„)„ converges to 

n h 

0 . 

2(p — 1) 

2. For p > 2 and q < -, we have 

p 

(a) {an)n converges to 0. 

(b) lim 

n —>+00 On 1 “t“ E 


(c) lim 


u. 


n+l 


n->-+oo U- 


= 1 + r > 1. 


Proof. — First of all, we look for some useful relations between a„ and a„+i. We recall ([^ 
for p > 2 and 1 
that (|^ implies 


2 (^p — 11 2 ^ 

for p > 2 and 1 < g < - < -We use the same calculations as (|12|) we obtain 


p 


p + l 


^ A„(l + r„(«” + (1 + 2A„)(1 + 


'^m-1 — 


1 + 3A, 


(14) 


Using (|5, we get 


u. 


®n+l — 


U. 


n+l 
m—1 
n+l 


1 + 2A, 


2Xnuit\ + (1 + Tniu:f,y-^)u:f, 


yn+1 


1 + 2 A^ 


2 An + 


(1 + r„K)»-i)'C ■ 


(15) 


u. 


n+l 

m—1 


By substituting (14) into (15) we get: 


fln+i < (1 + 2A„) < 2Xn + 


(l + 3A„)(l + r„«)P-^)n" 


< 


< 


A„(l + rn(n;^)P yu^ + {l + 2Xn)il + Tn{ul^_i)P 
A„(l + 2A.)(1 + + (1 + 2Xny{l + 

(1 + 3An + 2X‘fi){l + + 2An(l + 2An)(l + Tn{u^_fi)'P~yu'tf_i 

A.(l + + (1 + 2A.)(1 + 

(1 + A„)(l + r„(n" + 2A„(1 + 

An(l + rn«)P-^) + (1 + 2A„)(1 + Tn{u^n^-lY-yan 
(1 + An)(l + En(M5^)P“^) + 2A„(1 + Tn{vfn_f)P~yan 
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r 


But we have 7 ^— 7 , then 

'^n(l + t) + (1 + 2A„)(1 + r(a„)^ 


®n+l ^ 


(1 + A,i)(l + t) + 2An(l + T{an)^ ^)Ojr, 


And hnally we get 


®n+i -^n(l +'7‘)(an) ^ + (1 + 2An)(1 + r(a„)^ 

Un (1 + A,i)(l + r) + 2A„(1 + T{an)^~^)o,n 
In the other hand, using ([^ and Q we get 


n+l ^ + {I + Tn{u^_iy ^X_i r„ ^ ^ ^ 

1 + 2A„ h^(l + 2A„)^“'"^ 


n NQ-l n+l 


By using (|^, we have 


7/”+^ > 
“m-l — 


+ A„(l + + (1 + 2A„)(1 + 

(1 + 2An)^ 

2Anrn«)'?-^C+_\ rn(l + 


/l^(l + 2An)2 


h^(l + 2A„)2 


which implies 


2A2 , 2Anrn«)'?-' 


+ 


u, 


n+l 


> 


(l + 2An)2 ' h^(l + 2An)V 
An(l + rn«r^X + (1 + 2An)(l + rn«_,r ^)7i:;_, 

(1 + 2An)^ 

r^(i + rr.(uiy-y{uiy 


and then 


..n+l > 

«m-l — 


h^(l + 2An)2 

Anh^(l + rn«F-^X + h^(l + 2An)(l + TnK.i)^-^)71- 
h^(l + 2An)2 - 2hlXl + 2Anrn«)^-i 
n(l + r„«)>’-‘)(«”)» 


hn(l + 2An)^ — 2hnA^ + 2AnTn('U5^)'^ ^ 
Using (|^ and (18), we get 


®n+l 

yn+1 

^m-1 


— (1 + 2An) S 2An + 


(1 + 

«S-\ 


-1 


(16) 


(17) 


(18) 


> 


(1 + T„«)''-‘)< (\„hi - r„«)«-‘) + ft’(l + 2A„)(1 + T„(u” 


2\J&{1 + T„('U" +/!?(! + 2A„)(1 + r„('u;i)P-i)'ii;i 
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Then we can deduce that 

a„+i ^ (1 + T„«r ‘)(A„ - + (1 + 2A,.)(1 + 

On (1 + + 2An) + 2An(l + Tn{ul^_i)P~^)an 

Finally we get 

«n+i ^ (l + r)(A.-rh-nO^->-^ + (l + 2A.)(l + r(aOP-^) ^ 
o-n (1 + 'r)(l + 2Ajj) + 2A„(1 + r(a„)^’“^)a„ 


( 19 ) 


Next, we prove that the sequence (a„)„ 
need to show that < 1. 


But 


Cln 


converges to 0. 


To prove convergence, we only 


®n+i ^ A,i(l+ r)(an) ^ + (1 + 2An)(l + r(an)^ 

O-n (1 + A,i)(l + r) + 2A„(1 + T{an)^~^)o,n 


Let 


a4 (1 + r)A„ + (1 + 2A„)(1 + r^ttnY ^)cin ~ (1 + -^n)(l + 'T)(^n ~ 2A„(1 + T^ttnY ^)®n- 

We shall prove that A < 0. 

(1) First of all, we can see that, for p = 2 and g = 1 

A = (1 + t)A„ + (1 + 2Aji)(l + TaYjcin ~ (1 + -^n)(l + T)o-n ~ 2A„(1 + raYjoY 

= (1 + T)\n + (1 + Tan)(ln + 2An(l + Tan)(ln ~ (1 + ~ An(l + Yttn — 2A,i(l + TaYo-n 

= An (l + r + 2an(l + raY - (1 + Y^-n - 2a^(l + raY) + ral - ran 
((f T ^n) T 2(Xn(l T '^^YY ^n)) T f) 

= (1 - an)An(l + r + 2an(l + raY) + ™n(an - 1) 


Using 


we get 


I 0 < On < 1, 

ttn < U^_i, 

^ = ^nUY, 

^ hn 

A (1 (Xn)An(l T T T 2(Xn(l T raY ra^ii^a^i 1) 

< An(l — dYY + + 2(1 +!"))+ rn'lYn_i{an — 1) 

< (1-a„)rn(3h“^(l + r) 
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Using the condition: sup > 3h ^(1 + r),we can see that A < 0, so that 0 < a„+i < 

n 

On < 1, which implies that lim an = a exists and satishes 0 < a < 1. 

n^+oo 
2p-2 

(2) For p > 2 and q < -, we can see that 


p 


(1 “1“ -^n)®n ^ 0, 


if not, then, 


A„ r22-9 2-2p+pq 

+ >1 + Xn, 

Q-n '^m-1 


which is a contradiction because of o < 

^ p 

Let now. 


and 


Then it is clear that 
Ai > anA2. 


j4o — 


A 1 + T 
Ai = -T > 1 

^ 1 I p—1 

1 + ran 

‘2\n(^n ~ (1 + 2An) 


Xn (1 T Xn)ar 


< 1 . 


1 T T ^ (^ni^XnO^n (1 T 2A^)) 


1 T ToXn Xn (1 T XnXjOin 

(1 + T){Xn — (1 + Xn)cin) < fln(l + ^)(2Anan — (1 + 2An)). 

(1 + t)A„ + (1 + 2A„)a„(l + ra^ ^) — (1 + 'r)(l + A„a„) — 2A„a^(l + ra^ < 0. 
^ AI < 0. 

So 0 < ttn+i < a„ < 1. 

2ip — 1) 

We shall prove now that a = 0 for all p > 1 and 1 < q < -. By reduction to 

p 


absurdity we suppose that 0 < a < 1. Letting n —)■ cx) in (16) we obtain 

1 + raP-^ 

a < - a < a 

1 + r 

which is a contradiction. This proves that a = 0. 

AT 1 1- ®n+l Ip 1 2(p— 1) 

Next we prove that iim -= -, tor p > 2 and q < 

n—^+oc (Xn 1 T 

By means of we get 


p 


O‘n+1 ^ A„(l + r)(aTi) ^ + (1 + 2 Ati)( 1 + T(an)^ 

(1 + A,i)(l + r) + 2 Ati(1 + T{an)P~^)cin 


< 


( 20 ) 


but 


-2p+pq + 2 


At,( 1+ r)(a,,) =ci{ul^_^) (m”) 2 -? , 
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, r(l + r) 

where ci =-j-. 

And for q < — -we have — P + P? + — ^ then we obtain 


P 


2-q 


O'n+l , 1 

hm - < 


n— >+00 1 A E 


(21) 


In the other hand, using (19) we get 


^ (1 + r)(A„ - + (1 + 2A0(1 + r(a„)^’-i) 


> 


(1 + 'r)(l + 2A„) + 2A„(1 + r(a„)P 


but 




where Ci, C 2 G M. 

A j f 2(p - 1) 

And for q < -we have 


P 


(A„ - P) ^ 0 as n ^ + 00 , 


then we obtain 


O'n+l ^ 1 

hm -> 


—1-+0O Qj^ 1 A E 


( 22 ) 


Finally from (21) and (22) we deduce that 

lini ^ = ^<1. 

n^+oo dy^ 1 T 

To finish the proof of Lemma 2.3 we shall prove that for all p > 2 and q < 


2(p - 1) 


we 


have lim 


M, 


n +1 


p 


n— >■+00 nt 


= 1 A r. 


From ([^, we know that 

(1 + 2A„)«"+' - 2A„<y = (1 + r„(«” )>’-')«" , 

which implies 


7/”+^ 7; 

l + 2An-2A„^ = (lAr.(77::,r-)^ 


n 

m 

n+1 ' 
"m 


Then 


So 


1 A 2A„ — 2A„a„+i — (1 A t) 


K 


n+1 


U 


1 = (1 A r) lim 


K 


n—>-+oo U 


n+1 


U. 


n+1 


This implies lim —— = 1 A t > 1. This achieve the proof of Lemma 2.3. 


□ 


n-s>+oo u. 
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Now we can finish the proof of Theorem 2.2. We have showed that 


«m -2 < ( «m -2 + 11 


k=0 / k=0 


where 


An — 


(1 + r„<_ 2 ) 


1 + Aj- 


B„ = 


An(l + ^nU^)u^ + A„(l + 2An)(l + TnUl^_^)u^_^ 


(1 + A„)(l + 3A„) 

Using that >> 1, we can see that for p = 2 and g = 1 we have 


Then 


Tn = — and hn = h. 
uZ 


An < l+TnUl_, = l+T^ < 1 + = 1 + ra„, 


Bn = 


^ An(l + ^nUDu^ + An(l + 2A^)(1 + TnU^_^)u^_^ 
(1 + An)(l + 3An) 

^ A^(l + TnU^)u^ + A„(l + 2A„)(1 + TnU'^_i)u'^_i 

< A^(l + r)-u” + A„(l + 2An)(l + t„-u” 

^.n ^.n 

< c\l + r)^ + c(l + 2c)(l + r)^ 


with c = —. 
But we have 


< c(l + r) (c + (1 + 2c)) 

..n 

<c(l + r)(l + 3c)^, 


lim < X and > 0, 

n-5>+oo an 


0 < ttn < + 00 . 


then 
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In the other hand, for all c > 0, we have ^ ca„ < +oo, then 

n>0 

i<n(i + Cttn) < +00. 
n>0 


We deduce from this that 

0 < ^ < c(l + r) (1 + 3c) ^ ttn < 

n>0 n>0 


+ 00, 


and 


which implies that 


1 < JJ < JJ (1 + ran) < + 00 , 

n>0 n>0 


lim _2 < + 00 . 

n^+oo 

Now we will prove the second result of Theorem 2.2, that is: 


In we put j 


= m 


2(p — 1) 

If p > 2 and q < - - then lim m" _i < +oo. 

P n^+oo ^ ^ 

1 and we consider the quantity 


7;’^+! _ 7/"- 


X CA-xiAJ. vv^ ^v_/iioiLX\:::x aj^u-cixxlix ej 

- _1 < K{u^-2 - + Tn{u^_J 

= Gn + Hn, 


where 


G„ = kk:_\- 2 u';:_\ + kP) 

= c(«” - 2«”t\ +«”+‘) 

_ I n \-p+2+ ^),T^^ / '^rn-2 r, , i\ 

~ , n I n+l ^n+l 2a,i+i + 1) 

“m '^m-l 

2a-2 7 /"’+l 

“m “m-l 

> 0 , 


with c := —^ 
2 ^ 


and 


Hn = 


Therefore, using Lemma 2.3, we get 


= Tniu^.iY = rul^ittnY > 0 . 




^n+l . s-p+2+^ ^ , r /2(p-l) 

lim —— = (1 + r) 5-2 < 1 tor g < -, 

n^+oo Gji P 
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which implies 

Gn < + 00 . 

n>0 

Also 

lim = (1 + < 1 for p > 2 , 

n-S>+oo Hr, 


which implies 

< + 00 . 

n>0 


Hence we get the boundedness of 

“m- 

^ from: 

0 < <-i 

= 

n 

(“m-l - 



k=l 


< 

(Cfc-i+ Hk-i) + u\ 



k=l 



+ 00 


< 

(Cfc + Hk) + U^rn-l 



k=0 


< 

+ CX). 


Thus we have completed the proof of Theorem 2.2. □ 


3. Convergence 

In this section we prove the convergence of the numerical solution given by ([^, to the 
nodal values of the solution u of Q on each hxed interval time [0, T], T < T* as far as the 
smoothness of u is guaranteed. 

Lemma 3.1. — Let u be the classical solution of ([^ and he the numerical solution 
of Let T be an arbitrary number such that 0 < T < T*. Then there exist positive 
constants Co, Ci, depending only on T and uq, such that 

2(p — 1) 

(A) For p > 2 and q < - 

p 

max \uf — u{xj,t"')\ < Coh^~'^ 

l<j<m-2 ' ' 

holds so far as tn < T. 

(B) For p > 1 and q = 1 

max \u^ — u(x^,t"')\ < Cih"^ 

holds so far as tn < T. 

Before studying local convergence, we prove the consistency of the scheme. 
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3.1. Consistency. — For all 1 < j < Nn, we define 

u{Xj,t^^^) —u{Xj,t^) u{Xj+i,t"'^^) — 2u{Xj,t'^^^) + u{Xj_i,t"'^^) 


= 


T„, 


hi 


{u{x,,nY+ 


M(a:j+i,C) - u{xj-iA'^) 

q-l 

M(xj+i,C+^) — n(xj i, 

2hn 


2hn 


We use Taylor formula, we obtain 

u{Xj, — u{Xj, C) Tn d'^U 


du 

di 




du, 


Tn 2 dt‘^ 

u{xj+i, f”) — u{xj_i,tY d^u 


2K 


3 dx'^ 


(a;j,r + Tn9i). 
{Xj + hn02, C) 


(23) 


hl^ d^u 
3 dx^ 


{Xj - hnOz.tY- 


(24) 


d'^u 

dx'^ 


d'^u 

dx'^ 


{xj,r) = 


{xj,r) = 


u{xj+i, f”) — 2u{xj, C) + u{xj-i, C) hl^ d‘^u 


hi d% 


hi 


{Xj - hn05,tY- 


24 dx"^ 


Xj + hnOA, t") 


24 dx^ 

^YY(x- t'^+Y - T ^ ^ (x- + T Op) 

u{xj+i, t”+^) — 2u{xj, C+^) + u{xj_i, C+^) hi d^u 


+ + hn07,t'^^Y 


hi 


2Adx^ 


2Adx^^ ^ 

where 0 < < 1 for i = 1 ,..., i 

We define 


” -{Xj - hnOg.h^^Y + TnT-TTTjiXj.h^ + TnOp). 



du, 

<? 

F = 


— 


U 


dtdx"^ 


{Xj+i,tY - u{Xj i,tY 


(25) 


2hn 

We use the mean value theorem, the monotony and the symmetry of the exact solution 

Qll 

proved in Theorem 2.3 and Theorem 2.4 in [5], then there exists A between —{xj, f”) and 

such that 


2hr. 


du, 

9 


— 


u 


(a;j+i,r) - u{xj-i,f") 


2hn 

.n^ u{Xj+iA")-u{Xj_AA") 

q\A\ t ) 

= q\Ar^o(hl), 
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with 


Since 




< 


^n^ U(xj+ur) -U(xj-i,r) 

dx^^’ ^ 2K 


< o{hl). 


du 


dx 


So we can write that 

du 
dx 


is bounded before blow up by [T], then we can deduce that A is bounded too 




u 


u 


u 


u 


u 


2hn 
(xj+i,r) - 

2hn 

(a;j+i,r) - M(xj_i,r) 
2hn 

{xj+i.f') - M(xj-i,r) 

2hn 

(Xj+I,^) - u{Xj_i,r) 

2hr,. 


+ o{hl) 


+ o{hl) (26) 

2hn 

^{x„n+o{hi) +o{hi) 

du ^ 


9-1 


dx 


{Xj,e^^) +0{Tn)+0{hl). 


Then 



9 


9-1 



2hn 


2hn 



+o{Tn) + o(h^)- 




(27) 


We replace (23), (25) and (27) in e" we obtain 




d^u 


y + '^nOl) - ^{Xj, n - r + rj^) 


h‘^ h‘^ 

- hn9e,r^^) - Hx„nr 

r\ Q 

+ o{Tn) + o{hl). 


+ 


du, 


If we put 


Ri = 


max 

X.t 


1 d‘^u. , d^u . , 


and B.o = 


d^u 


we can deduce that 

max e” < CiTn + C' 2 h^, 

with Cir„ = RiTn + o(r„) and C' 2 h^ = -^ 2 ^^ + c'(^n)- 
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3.2. Local convergence. — Let e” = m” — u{xj,t"') for j = 1, ...,m — 2. 
(A): Using (23), (25) and 26 we get 

u{xj, — u{xj, t"-) — 2u{xj, + u{xj_i,t"'~^^) 


Tn 


hi 


-Hxj,nr 


+ 

Tn d‘^U 


U 


(xj+i,r) - u{xj-i,t'^) 


2hr. 


d^u 


2 df^ 


{Xj,r + OiTn) - t" + 6 * 6 ^ 71 ) 

h‘^ f f)^l! 

-d 7 <1 + OrhnX^^) + - eshn,r^^) \ + o{K). 


Let 


r” •= 

3 ■ 


24 \ dx‘^ 

Tn d‘^U 


2 df^ 
hi fd^u 


d^u 

{Xj, r + 9iTn) + Xn Q^Q^2 + ^e^n) 


d% 


1 ^ ~ )} + o{hJ 


dx* 


Then 


u{Xj, — u{Xj, t”) u{Xj+i, — 2u{Xj, + u{Xj-i, 


hi 


+ 


u{xj+i,r) - u{xj_i,r) 


2hr. 


= —T\ 


-Hxj,nr 

(28) 


Using (|^ , we have 


{Uj) -\- 


hl 


{2hny 




= (29) 


From (28) and (29), e" satishes 


n+l _ n n+1 _ 9 „n+l , n+1 

i ^3 U+1 ^ U-1 


hi 


-{{u-r-u{x,Xr) 


{2hny 


« . Aq-Mu^+I - _ 


m,, 1 - u 


i+1 U-ii lU+i U-i 


M 


(a:j+i,r) - M(a:j i,r) 


2 h, 


_ 

i • 


By the mean-value Theorem, for f{X) = X^, we get 

{u]y - {u{x,,ny = f{uy)-f{u{x,,n 

= f{u{x,x) + e,eyey 
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for some 6g G [0,1]. Then we obtain 

^n+l 

"3 '-'3 "i+1 “V 


e?+^ - e? e?+^ - 2 e^+^ + 


Tn hi 

= f'iu{xj,r) + ege])e] 


{2hnY 


\u 


i+l '“i-l 


+ 


u{Xj+i,t^) - u{Xj-i,t^) 


2hn 


+ r;. 


l“i+i 


u 


n+1 1 

'i-i I 


Using (26) we get 


n+1 _ n n+1 _ o^n+l , n+1 

S' S+1 ^ S-1 


hi 


n 
3 


= /(«(x,,r) + «)e 


with = r” + o{hl). 
Let: 


( 2 h„)^ 


1 'n 


n |<J-1 l^^+l _ r,,n+l 


'j +1 “j-ll l“i+l “j 


wrii + 


-(x- t 


n^ I „n 


E^= max |e"|, 




U = max |■u(x, t)|, 


X,t 


V = max 

X^t 


du, 


W = - max 

3 a:,i 


d^u 


dx^ 


x,t) 


K = f{U + l), 


and 


D ''n 

K = — max 

2 x,t 


d‘^u 


x,t) 


+ max 

X± 




dx'^ 


x,t) 


H-max 

12 x,t 


d% 

dx^ 


x,t) 


+ o(l) + o(A„). 


But from (24) we have 


du, , 


“i+i “i-i 


2hr. 




“i+1 “i-i 


2hr. 


3 dx^ 


2 hr, 




h^ d^u 


2 hr, 


rpn 

< J- + hlW. 


3 


(X, + hhrr,r) - y - ehn,n 


(30) 


Then by (30) and the mean value theorem, for g{X) = |X|'^, we get 

q 


du , 

q 


~ 


“i+i “i-i 


/ f du 


2 h„ 


- ''S ( + —+ '<ir — + KW 



E” 


(31) 
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In the other hand, for 1 < j < m — 2 we have, 


“i+1 

q 

“i+i “i-i 

q-1 

.,n+l _ n+1 

“i+1 ^i-1 


2hn 


2hn 


2hn 



< 


< 





1-^ ( r,,'^ — 

j + 1 “i-1 


^j + 1 


D -1 


2 hr, 


2hr 


9-1 _ pfi pn+1 _ n+1 

j + 1 ^j-1 ^j+1 ^j-1 


2hn 

q-1 / pU _ pU 
j+1 ^j-l 

2hr,. 


2hr 


+ o{Tn) + o{hl) 


+ 


“i+1 “i-1 


2hr 


q-1 / n+1 _ „n+l' 

^i+1 ^i-1 


2hr 


+ {o{rn) + o{hl)) 


^j+1 - ^j-1 


2hr 


q-1 


(32) 


Then from (32) and (31) we get 

0^+1 _ 9 „n+l I n+1 


p^+1 _ pn p- 


hi 


< f{u{xj, r) + Oge^)e^ + r” + (o(r„) + o{hl)) 


^j+i - ^j-i 


2hr 


q-1 


J I 



hr. 


“i +1 “i-i 


2hr 


q-1 /pU 


"i +1 t -1 

2hr. 


Let M := ||?7”||oo = Finally we obtain 

^n+l _ 


q-1 


< KE'' + hlR+(o(+) + o(hl))[- 


'M' 


9-1 


' (^ 1 / + 0 


jpn 



+ + (o(An) + o(l)) J +Wqg 

rpn 

= JT (KX + 13' {+V + 0 (£” + hlw ))) 

+ I’jr'' {ht'R + (o(A„) + o(l)) M®-' + {h„V + 9 (B“ + AjW'))) 
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Let 


Then, 


B = hlK + qg'{hnV + e {E^ + hlW)) 

C = hl-^R+{o{\n) + o{l))M‘i-^+ Wqg'{KV+ e{E^+ hlW)) 


< (^l + Tn^^E^ + Tr.hl-'iC 

< {l + TnNB)E^ + Tnhl-‘^C 


< exp{NBT)hl-'^CT. 


With N is constant snch that 


»7(g-l) 

1 I ][f 2-q 

For tn<T and we have: ^ := iV, 


hi 2 ^ 


which is bonnded by Theorem 2.2. Then we get 


3-g 


max \u'l — uixj,E)\ < Co{T)h 
Now, we will prove the last part of the lemma. 

(B): We do the same thing for p > 1 and g = 1, we get for j = 1,..., m — 1 


„n+l _ 
^3 


.. n+1 

S' S +1 


2^+1 ^ gn+l 


"'n 


= f'{u{xjX) + 09e^)e^- 


S +1 S -1 

2h„ 


+ 


u{xj+i,E) - u{xj_i,r) 


2K 


+ r" 


gn+l _ gTi+l 

= f{u{x„E) + ege])e]+ + r ■. 


And then 


2 h^ 


^n+l _ 


< KE^ + hlR. 


And hnally we obtain 


+ TnhlR. 

< exp(A:T)h^i?T. 


max In” — n(x,', S) I < (FifT)/;,^. 


4. Approximation of the blowing up time 

In this section, we give an idea abont the nnmerical blow-np time. First of all we recall 
a resnlt of Sonplet and Weissler |S] 
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Theorem 4-1- — Let -0 G (s large enough), with ‘ip > 0 and '0 7 ^ 0. 

1. There exists some Aq = Ao('0) > 0 such that for all A > Aq, the solution of Q with 
initial data (p = \fj blows up in finite time in norm. 

2. There is some C > 0 such that 


T*(AV’) < 


C 


A —)• 00 . 


3. 


T*{\ij) > 


1 


We define now 


T* ■= T 

num * / j ' '> 

n>0 


(33) 


and call it the numerical blow-up time. In [5], we have proved that 


ul > 


l + r 


-q , „ ^ - 2 p+ 9 (i+p) 
. 1 -h r22-9 2-5 


n 


U 


0 

m* 


which implies that 


< 




l + T 




(34) 


a + r22-9(M^) 


-2p+g(l+p) 

2-9 


Using (33) and (34) we get 


rri* 

^ num 


n>0 



/ 


< 


T 



E 

n>0 


\ 


1 


_l+r_ 

-9 -2p+9(l+p) 

l+r22=9 (u^) 2=9 


\ 


r 



E 

n>0 


1 -|- r22-9 




-2p+9(l+p) 

2-9 


l + T 



n 


T 



l + r22-9 (n)),) 


-2p+9(l+p) 

2-9 


p-1 


_ 


(35) 
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5. Numerical simulations 


In this section, we present some numerical simulations that illustrate our results. In 
hgure 1, we take p = 4 > 2 and g = 1.3 < one can see that the solution is bounded 

in Xm-i- Then we take p = 2 and g = 1, it is clear from hgure 2 that the solution blows 
up in Xm-i, and from hgure 3, we can see that the solution is bounded in Xm- 2 - 
Concerning the approximation of the blowing up time, if we take the initial data Uo(x) = 
Asin(|(a; + 1)), with A > 0 then ||mo|Ioo ~ Theoretically we know that 


T* > 


(p-1) \\uo\\ 


P-1 * 

OO 


Let g(A) = 


—- and p = 3. In the next table, and for some values of A we can 

D—1 ^ ’ 


{p — 1 )Ap“ 

see that > g(A) which is compatible with the theoretical result, this is illustrated in 
hgure 4. Also, using (35) and for A = 10^ we have 


T* ~ 
num 


5.067.10“^ < T** = 5.075.10 


-7 


A 

10 

10^ 

103 

10^ 

103 

gW 

5.10-3 

5.10-3 

5.10-^ 

5.10-9 

5.10-1^ 

rji^ 

num 

5.177.10-3 

5.068.10-3 

5.067.10-^ 

5.075.10-9 

5.058.10-1^ 


Table 1. Comparison of the function g{X) with the numerical blow up time T*, 



Figure 1. Evolution of the 
numerical solution at Xm and 
Xm-i for p = 4 and q = 1.3 



Figure 2. Evolution of the 
numerical solution at Xm and 
Xm-i for p = 2 and g = 1. 
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Figure 3. Evolution of the 
numerical solution at Xm and 
Xm -2 for p = 2 and q = I- 



Figure 4. Graphics of g{X) 
and approximation of the nu¬ 
merical blow-up time for p = 3. 


6. Conclusion 


We have showed that when p = 2 and q = 1, the hnite difference solntion blows np at 

2(p — 1) 

more than one point and that when p > 2 and q < -, the only nnmerical blow np 

p 

point is the mid-point x = 0. This is an interesting phenomena in view of the fact that the 
solntion of the corresponding PDE blows np only at one point x = 0 for any p > 1 and 

2p 2(p — 1) 2p 

1 < q < -. Remark that for 1 < p < 2 and - < O' < -, we have no idea 


p + 1 

abont the bonndedness of n^-i '^m- 2 - 


p 


p + 1' 



Blow-up of the solution at 
Bounded solution at Xru—\ 


Figure 5. Graphics of the asymptotic behaviours of the solution near the 
blowing-up point. 
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